Nuclear applications of inverse scattering, 
present . . . and future? 



R.S. Mackintosh 
Department of Physics and Astronomy 
The Open University 
MKT 6AA, UK 



Abstract 

There now exists a practical method (IP) for the routine inversion of S- 
matrix elements to produce the corresponding potential 0. It can be applied 
to spin-i spin-1 projectiles. We survey the ways that the method can be 
applied in nuclear physics, by inverting Sij derived from theory or from ex- 
periment. The IP method can be extended to invert Sij{E) over a range of 
energies to produce a potential V{r,E) + l-aVis{r, E), yields parity-dependent 
potentials between pairs of light nuclei ,2 and can be convoluted with a di- 
rect search on the S-matrix to produce 'direct data — > V inversion'. The 
last is an economical alternative form of optical model search to fit many 
observables (e.g. for polarized deuterons) for many energies, producing an 
energy-dependent potential with many parameters (e.g. Tr for deuterons) [3]. 



1 Introduction 



It is very easy to derive a cross section or an S'-matrix from a potential; the re- 
verse is much harder, but can now routinely be achieved for a wide range of cases. 
Here we introduce the IP (iterative-perturbative) inverse scattering method, with 
the emphasis on its range of applicability, illustrated by a few diverse successful 
applications. We hope to inspire applications that we have not thought of. 

Alternative methods exist for some tasks to which IP inversion can be applied (e.g. 
conventional OM searches to obtain potentials from scattering data, and weighted 
trivially equivalent potentials (TELPs) to derive DPPs) but IP inversion not only 
has advantages in these cases but also gives reliable results where no other techniques 
are available. 



2 Inverse scattering: problems and solutions 

'Inverse scattering' usually refers to the derivation of the potential corresponding 
to given S'-matrix elements or phase shifts; it may be contrasted with the trivial 
forward case. This is Case 1 of Table ^ We shall say little about Case 2 of the 
table, but we will mention Case 3, direct inversion from observables to potential, 
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since the IP method can readily be convoluted with Case 2 inversion to yield a very 
powerful Case 3 alternative to conventional optical model searches. 



Case 1 


forward 
inverse 


V{r) ^ Si 
Si ^ V{r) 


EASY 
MUCH HARDER 


Case 2 


forward 
inverse 


Si ^ a{e) 
<j{e) ^ Si 


TRIVIAL 
OFTEN VERY HARD 


Case 3 


forward 
inverse 


V{r) ^ a{9) 
a{e) ^ Vir) 


EASY 
VERY HARD 



Table 1: Traditional characterization of forward and inverse scattering cases. The 
a symbolizes all observables (analyzing powers, etc); Si includes Sij etc; V includes 
spin-orbit and tensor terms. 



2.1 Traditional methods of Si — > V{r) inversion 
Formal inversion methods apply to two classes: 

1. Fixed-/ inversion. (Gel'fand and Levitan; Marchenko) Si — > V{r) for Si for a 
single / for all energies, yielding a local potential. 

2. Fixed energy inversion. (Newton-Sabatier (NS); Miinchow and Scheid (MS)) 
Si — > V{r) given Si for all I at a single energy. Related procedures have been 
developed by Lipperheide and Fiedeldey. 

There also exist semi-classical methods for fixed-energy inversion based on WKB 
and Glauber (eikonal) methods. Practical versions of NS, MS, can handle a finite 
range of /. The review [T] gives comprehensive references to the inversion techniques. 

Disadvantages of traditional methods: 

(1) The NS method requires highly precise Si, and there is a tendency to instability. 

(2) The Marchenko method requires large energy ranges, but nuclear potentials are 
energy-dependent (but it has been used for nucleon-nucleon scattering). 

(3) They are not adaptable to cases with small ranges of / (NS). 

(3) They mostly apply to spin 0; NS can handle spin 1/2. 

(4) They are not readily generalizable. 

Nevertheless, there have been a number of applications which have yielded real 
physical insights, see Ref. p. 



2.2 Information on nuclear interactions from inversion 

Inversion can be applied in three distinct general ways: 
I. Inversion of Si, Sij or S^, obtained from theory: 

1. Derive dynamic polarization potentials (DPPs) arising from (i) Inelastic scat- 
tering, (ii) Breakup processes, (iii) Reaction channels, etc. 

2. Derive potential from RGM and similar 5-niatrix elements. 

3. Determine local potentials 5-matrix equivalent to non-local potentials. 
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4. Obtain a potential representation of impulse approximation S-niatrix, or the 
5-matrix from Glauber model and other impact parameter models. 

II. Inversion of Si, Sij or S^, from analysis of experiment: 

1. (When there are few partial waves.) Invert Si from parameterized _R- matrix 
or effective range fits at low energies. (Requires 'mixed case' or 'energy de- 
pendent' inversion, see Section 3.3) 

2. (When there are many partial waves) High energy two-step phenomenology. 
(E.g. "Li, +^^C from 140 to 2400 MeV) 

III. Direct observable V{r) inversion. 5-matrix search can be convoluted 
with IP inversion yielding the 5- matrix is byproduct. Many energies can be treated 
simultaneously to give a multi-component V{E). 



3 The Iterative-Perturbative (IP) method 
3.1 The key idea 

The response of the elastic scattering S'-matrix to small changes is assumed to be 
linear (this is often surprisingly accurate): 

^^' = -5i: / iMr)rAV{r)dr. (1) 
n k Jo 

where ui(r) is normalized with u/(r) Ii{f) ~ SiOi{r), Ii and Oi are incoming and 
outgoing Coulomb wavcfunctions. Eq. ^is readily generalized for spin. 



3.2 An Outline of the IP method: 

Take a known 'starting reference potential', SRP, V{r) giving Si. With added term: 

V{r) ^ Vir) = Vir) + c^v,ir) (2) 

it gives Si + ASi. Fimctions Vi{r) belong to a suitable 'inversion basis'. 

The core of the method is the solution using SVD of the over-determined linear 
equations derived from Eqn. ^with ASi = S^'^^^^^ — Si and AV = X]c^^'^(^) to 
find amplitudes Ci such that V gives Si closer to S^^^^^^ . By iterating the linear 
equations, Si -I- AS"; converges to S**^"^^"*. There is often a natural starting potential; 
it can often be zero, or the 'bare potential' when establishing DPP contributions. 

The facility of the IP method to control the fitting is a key element in the method, 
especially in view of the innate ill-posedness (Ref. ^) of the inversion problem; 
it is possible always to demand smooth potentials. This and all aspects are fully 
discussed in Ref. PP with many references. 



3 



3.3 Generalizing from fixed-energy inversion 



IP Inversion can be generalized indefinitely as the following progression suggests: 

1. Fixed-energy inversion Si, 'all /, one E' inversion. However, at low energies 
the potential is generally under-determined, there being too few active partial waves 
to define the required potential. 

2. Mixed case (energy bite) inversion. The problem of under-determination 
at low energies can often be solved given Si{E) over a range of energies ('energy 
bite'). This is 'some /, some E\ Si{E) V{r), or 'mixed case', inversion. For a 
narrow energy bite, this is effectively includes dSi/dE as input information. 

3. Energy dependent inversion. Nuclear potentials, particularly the imaginary 
parts, vary with energy, but IP inversion can be extended to determine V{r, E) 
directly: 'some I, some E\ Si{E) V{r,E). To date, cases have been limited to 
the factored form fi{E)Vi(r), where i indicates real-central, imaginary-central, 
spin-orbit, Majorana, etc. 

4. Inversion to fit bound state and resonance energies. Energies of bound 
states can be included with Si as input information to determine V. 

5. Direct data to potential inversion. Example: d— *He, multi-energy. 

IP inversion can be applied to the case of identical bosons where only even partial 
waves are involved, e.g. ^^C -I- ^^C. It also applies to the case, very important with 
light nuclei, where the potential is parity dependent. Such a potential can either 
be represented as a sum of independent even-parity and odd-parity terms, or as a 
sum of Wigner (W) and Majorana (M) terms: 

Vw{r) + {-lfVM{r). 

It is always found that even-parity and odd-parity potentials have different radial 
forms which often implies a surface peaked Majorana term. Both RGM ^-matrices 
and experimental data imply that there is significant parity dependence even for 
nucleons on nuclei as heavy as ^®0. 



3.4 Spin cases that can be handled by IP inversion 

1. Spinless projectiles. Si — > Vir)- 

2. Spin 1/2 projectiles. 5*;^ V{r) + 1 ■ crVis(r). 

3. Spin one projectiles. Vector spin-orbit and Tr = ((s • f)^ — 2/3)VR(r) tensor 
potentials can be determined from non-diagonal 5/;, . This is coupled channel inver- 
sion. 

4. High channel spin. For cases like d -I- ^He, independent potentials for each 
possible channel spin have been determined. 

In every case, all spin-dependent components may have real and imaginary and 
Wigner and Majorana terms. These can all be expanded in different bases. 



3.5 How well does it work? 

Fig. 1 shows a test case [H] in which a Sf^, for deuterons on ^^Ni at 56 MeV and a 
known potential, including a tensor term, were inverted with an arbitrary SRP. IP 
inversion can be applied to noisy data and produce meaningful interactions because 
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the departure of the final potential from the 'starting potential' of the iterative 
method is under control, see Section 4.3. 




Figure 1: Deuterons on ^^Ni at 56 MeV. Solid lines: target (known) potential; dots: 
V found by inversion; dash-dot: inversion SRP (starting potentials, zero for the real 
spin-orbit and tensor terms.) 



4 Selected applications of IP inversion 



Light nuclei, parity dependence. Scattering between various pairs of light nu- 
clei have been studied by inverting S'-matrices from RGM calculations and from 
i?-matrix fits to experimental data, over a wide range of energies. These stud- 
ies reveal the importance of a parity-dependent (Majorana) component. Ref. ^ 
has references to parity-dependent potentials for various pairs of light nuclei and 
presents a case study of p+'^He, comparing potentials from empirical and theoreti- 
cal Sij. The contrasting Majorana potential for d-|-^He is discussed in Ref. |6j and 
that for p+^Re in Ref. E]. 

The dynamic polarization potential, DPP. It is well known that the coupling 
to breakup channels generates a repulsive DPP for projectiles such as ^Li and ^H. 
Inverting the elastic scattering S'-matrix from a coupled channel calculation, and 
subtracting the bare potential, gives a local-equivalent /-independent representation 
of the DPP. The form of the DPP depends on the L-transfcr in a systematic way |3 . 
Nucleus-nucleus interaction also receive large contributions from coupling to inelas- 
tic and (especially) reaction channels that cannot be represented by renormalizing 
folding model potential. Many cases are described in Ref. _1_, and the contribution 
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of breakup to the p+^He interaction is presented in Ref. [S], revealing the hmitations 
of folding models for halo nuclei. 

Potentials from empirical data. Elastic scattering potentials, including ^'^C 
_l_ i'2(j^ 16q _|_ 16q^ iij^j scattering and nucleon scattering, have been determined 
either by fitting S to data and then inverting, or by using the 'direct inversion' in 
which the S'-matrix search is convoluted with the S V inversion. The definitive 
phenomenology for has been carried out in this way revealing, inter 

alia, the necessity for a Majorana term. Direct inversion of multi-energy data for 
scattering of light nuclei provides an alternative means of establishing phase shifts 
Sij{E) that behave in a way that is consistent with a potential 1101 that varies 
smoothly with energy. Inversion of Sij from R-matrix fits also has this property. 

4.1 Pickup coupling effect in ^He(p, p) elastic scattering 

It known that pickup coupling {p ^ d ^ p for proton scattering) makes a significant 
contribution to the nucleon optical potential. This is one reason that precise fits for 
nucleon elastic scattering below 50 MeV for closed shell target nuclei have not been 
found with conventional optical model fitting (for non-closed shell nuclei it is easier 
to find parameters that fit the shallower diffraction minima). It is now possible 
to do full finite-range pickup calculations including non-orthogonality terms, and 
these have been carried out [2] for the p-*He system. The coupling has a large 
effect on the elastic scattering angular distributions. In Table 2 we quantify the 
pickup contribution in terms of volume integrals. The repulsive real DPP is quite 
large at the nuclear centre although the effect on the volume integral is modest. 
The radial form of the DPP could not be represented by renormalizing a folding 
model potential. The volume integrals reveal the importance of including the non- 
orthogonality correction. 





Jk 




Jl 


/ 2\l/2 


JSOR 


■^SOI 


OM 


704.14 


3.092 


55.37 


3.336 


26.60 


0.005 


CRC 


653.94 


2.938 


307.47 


4.138 


40.27 


1.25 


NONO 


571.28 


2.840 


252.62 


4.360 


33.15 


6.55 



Table 2: For 15.6 MeV protons on *He, volume integrals per nucleon pair/(MeV 
fm'^), and rms radii/fm of the bare potential (OM) and the potentials found by 
inversion for the complete CRC calculation and for that in which non-orthogonality 
term was omitted (NONO). 



4.2 The DPP due to breakup for ^He scattering from ^°^Pb. 

The breakup for this case ^1 E| was calculated using the same model for ^He as 
Ref. [Hj. However the DPP is now very different, having a long range attractive tail 
generated by the Coulomb dipole interaction. Fig. 2 shows the short range repul- 
sive/emissive DPP in the surface region. Local regions of emissiveness are a common 
feature of local potentials representing the highly non-local and Z-dependent dynam- 
ical polarization contributions; unitarity is not broken. The DPP is not well-defined 
for r < 10.5 fm. 
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Figure 2: For ^He incident on ^°*Pb, 27 to 66 MeV, the DPP at the nuclear surface. 



Fig. 3 shows the long range attractive and absorptive potentials generated by the 
dipole coupling; the real part is appreciable out to 60 fm. Both the real and imagi- 
nary DPPs strongly influence the elastic scattering differential cross section. 



4.3 The d-^He interaction derived from multi-energy data 

'Direct data — > V inversion' is an alternative ^U] to optical model fitting with 
parameterized forms for determining potentials from data, especially when there 
are many data and many parameters. This is the case for d-'^He scattering when 
there is a full set of polarization observables (including all 3 tensor analyzing powers) 
for many angles, all for many energies ranging from 4 to 13 MeV. Refs. [3] fitted 
1000 data points (five observables, a wide angular range and many energies) to 
produce a multi-component (Wigner and Majorana, central, spin-orbit and tensor) 
multi-energy potential (components were functions of energy) giving a reasonable 
representation of shape resonances. This would have been a formidable task for 
standard optical model codes since they would have had to include the coupled 
channel calculation for the tensor interaction within the search. 



5 Possible future applications 



1. Systematic CRC calculations followed by inversion of the resulting elastic Sij 
provides a method for establishing shell corrections to the nucleon OM potential. 

2. There exists much data for the elastic scattering of neutrons from ^^C. Inversion 
of phase shifts fitted to this data would yield a potential model fitting the non- 
resonant scattering, allowing an extrapolation to higher neutron energies. Compar- 
ing the potential derived in this way with theoretical and standard empirical OM 
potentials would support (or the opposite) the neutron scattering data as well as 
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15 20 25 30 



r(fm) 

Figure 3: For ^^He incident on ^o^Pb, 27 to 66 MeV, the DPP outside the nucleus, 
the theoretical models. 

3. The IP algorithm appears to be indefinitely generalizable, and has found wider 
applications than originally envisaged, so are there new applications to nuclear data 
evaluation? Any suggestions? 
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